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Abstract. The transport coefficients of the Anderson model require knowledge of 
both the temperature and frequency dependence of the single-particle spectral densities 
and consequently have proven difficult quantities to calculate. Here we show how these 
quantities can be calculated via an extension of Wilson's numerical renormalization 
group method. Accurate results are obtained in all parameter regimes and for the 
full range of temperatures of interest ranging from the high temperature perturbative 
regime T >> Tk , through the cross-over region T « Tr- , and into the low temperature 
strong coupling regime T « Tk- The Fermi liquid relations for the coefficient of 
the resistivity and the linear coefficient of the thermopower are satisfied to a high 
degree of accuracy. The techniques used here provide a new highly accurate approach 
to strongly correlated electrons in high dimensions. 
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1. Introduction 

In this Letter we present accurate new results for tlie transport coefficients of the 
Anderson model obtained by using the numerical renormalization group method [Tj. 
The Anderson model has been used extensively to interpret the properties of magnetic 
impurities in simple metals, certain aspects of heavy fermion behaviour [2] and the low- 
temperature transport through a quantum dot [3] . Recently this model has also become 
important for an understanding of the properties of correlated electrons on a lattice in 
high dimensions. The electron self-energy becomes local in infinite dimensions |H |5] and 
the problem reduces to that of an impurity [6]. In particular the solution of the d = oo 
Hubbard model has been shown to reduce to the solution of an Anderson impurity model 
with a hybridization function, A(a;), determined by a self-consistency requirement 
[ZllHlin]. The self-consistency condition is straightforward to handle and the difficult part 
of the problem has remained the calculation of the dynamics of the underlying impurity 
problem, and more specifically of the impurity single-particle spectral densities. This 
is just the problem which one encounters in attempting to calculate the transport 
coefficients of the Anderson model which depend sensitively on the temperature and 
frequency dependence of the single-particle spectral densities pd{uj,T). An accurate 
solution of this problem can therefore provide the solution of the d = oo Hubbard 
model. Our main interest, then, is the accurate solution of the Anderson model and in 
particular the calculation of its single-particle spectral densities. The Bethe Ansatz has 
been successful in providing exact results on the thermodynamic properties of several 
models of magnetic impurities including the N-fold degenerate Anderson model [TOlITT]. 
but the method cannot be used to extract dynamic properties. The latter have been 
obtained by applying approximate methods. For N = 2 perturbation theory in the 
local Coulomb repulsion, U, has been used [12] whilst for larger N (typically > 6) the 
non-crossing approximation (NCA) has been applied to calculate the thermodynamic, 
transport and dynamic properties [131 El- The former eventually becomes unreliable 
for large enough values of U (typically for U/nA > 2.5 where A is the resonance level 
width) and the latter fails to satisfy some exact Fermi liquid relations at low energies 
[15] . Recently, Quantum Monte Carlo methods in combination with the maximum 
entropy principle have also been used to extract dynamic properties [lEl [ITJ [18] . These 
have so far been restricted to the symmetric model and to fZ/vrA < 3. They become 
more inaccurate and time consuming for larger values of U and lower temperatures. 

1.1. Model and method 

A comprehensive and highly accurate approach to the calculation of dynamic properties 
of magnetic impurity models has recently been developed, based on an extension of 
Wilson's numerical renormalization group approach, which overcomes many of the 
limitations and difficulties of the other approximate schemes. This has led to reliable 
results valid in all regimes and at both zero [HI [201 [211 [221 [231 [2l] and finite [25] 
temperature, and which satisfy all the sum rules and Fermi liquid relations of the 
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Anderson model [2U]. With such results for the dynamic properties it becomes possible 
to calculate the transport coefficients of the Anderson model over the full temperature 
range of interest. We report here results for the resistivity, and thermopower based on 
these numerical renormalization group calculations. Results for the Hall coefficient and 
thermal conductivity have also been calculated and will be presented in a forthcoming 
publication \26\ . 

Our starting point is the non-degenerate (A^ = 2) Anderson model: 

-f^ = ^ (^dC^^Cda + U n-d^Udl + '^{Vkacl^Cda + H.C.) + ^ efcc|.^Cfc^ (1) 
(T kcr kcr 

This can be re-cast in the form of a linear chain model which may then be iteratively 
diagonalized. The Lanczos algorithm is applied to He = ^kcl^Cka with starting 
vector CQ^\vac >= yj^k^^cl^vac > where V = ^/^^jWF so that the conduction 
electron part of H is reduced to a semi-infinite linear chain with site coupled to the 
impurity via a hybridization term of strength V [H [12] : 

-f^ = ^ edcl^Cda + U ridtUdi + V ^(^^q^ + H.c.) 

a a 

oo oo 

+ ^ ] ^nf^na^na + ^ ] \i{c'n^\fj'^ncj + ^na^n^Xcj) (2) 
n=0,a- n=0,cr 

The parameters {e„,A„},n = 0,1,... reflect the form of the hybridization function 
A(cj) = 7r^^|Vfcp5(a; — e^). In the present case of the magnetic impurity problem 
we are primarily interested in low energies where A(a;) may be taken to be a slowly 
varying function of lo which can be approximated by a constant A. In this case, and 
by replacing the continuum of conduction electron states in (— -D, U) by a discrete set 
±Z}A~",n = 0,1,..., A > 1, a discrete approximation to H can be obtained with 
parameters A„ ~ DAr'^/'^^n » 1 [1| (the e„ are zero for a half-ffiled symmetric band). 
For the impurity model corresponding to the d = oo Hubbard model, the frequency 
dependence of A(w) is important and leads to a different set of parameters A„. The 
model has to be solved repeatedly in this case to obtain a self-consistent A(c<;) using 
the self-consistency requirement in [TJ El E]. 

The discretized Anderson model with A„ ~ DA~"'^^,n » 1, discussed above, can 
be iteratively diagonahzed to obtain the many-body eigenvalues and eigenstates 
\p >n on successively lower energy scales Un ~ DA~"/^,n = 0,1,... by following 
the procedure in [1]. The procedure starts with diagonalizing the impurity part 
^0 = J2a ^dc\„Cda + Uud^Udi and then adding the coupling to the local conduction 
electron state V J^ai'^Ocr'^d.'^ + H.c). After this one adds successive energy shells 
Xlo- -^"('^n+io-Cncr + c'la'^n+ia) , IT- = 0,1,... and diagoualizes the resulting Hamiltonian. 
Since the number of states grows like 4" it is in practice only possible to retain the lowest 
500 or so states at each stage. This restricts the reliable range of excitations u to be such 
that < a; < Kun, where i^' ^ 10 for A ~ 3. The lower excitations are calculated 
more reliably in successive iterations, whilst information on the higher excitations is 
contained in previous iterations. The eigenvalues are used to calculate the partition 
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function at a decreasing sequence of temperatures T„ = Un/ks ~ -DA "/^,n = 0, 1, . . . 
and from this one extracts the thermodynamic properties [H |3T] . 

2. Results 

2.0.1. Spectral functions From the eigenvectors \p >„ we recursively evaluate the 
single-particle matrix elements iW^, =< > which are required for the n'th shell 

Green's function G^^i^uj.T) and spectral density p'^{u,T), 



1 ^M^re-'^^p+e 



Here Zn{(3) is the partition function for the n'th clusteil§| or energy shell. In evaluating 
p2{u!,T) the delta functions are replaced by Gaussians having widths an of order cOn 
appropriate to the cluster size. At T = the spectral density is evaluated at excitations 
w ~ using the n'th cluster. At finite temperature /c^T > 0, and when u ^ Un becomes 
of order ksT excited states not contained in the n'th cluster become important. In this 
case we use a smaller cluster (containing information about the higher excited states) to 
evaluate the spectral densities at energies u such that u < ksT. Results for the spectral 
densities pd{uj,T) are shown in Figure [Hand Figure [2] for various temperatures. 

The Friedel sum rule pd{uj = ep,T = 0) = -^^sin^ {1112^/2) and Shiba relation for 
the dynamic susceptibility [3l] are satisfied to within a few per cent in all parameter 
regimes [HI [201 135] . 

2.0.2. Transport coefficients The calculation of transport properties is straightfor- 
ward once accurate results for the temperature and frequency dependence of the spectral 
densities are obtained. The resistivity p(T) and thermopower S{T) are given in terms 
of the integrals 

Lmi = -C ^J^r\uj)u:^du (5) 



by 



PiT) = (6) 



01 



where the relaxation time t{u,T) is given in terms of the local d-electron Green's 
function by t{uj,T)^^ = 2n\V\'^ pd{ou,T), and e is the electronic charge. There are two 

§ The term cluster is suggested by the notation, ahhough the calculations described here are in k-space 

ED 
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exact Fermi liquid relations for the transport coefficients. The first is for the coefficient 
of the term in the low temperature resistivity |271 |2H]- For T « and in the 
Kondo regime we have 

P(T)=p(0)|l-c(^^y| (8) 

4 

where c = = 6.088, and Tk is defined by 



2U J 

Results for the resistivity are shown in Figure |3] for U / tt A = 4 and for several values 
of the local level position ranging from the Kondo regime (e^/A = — —4, —3, —2) to 
the mixed valent regime (cd/A = —1,0) and the empty orbital regime (e^/A = +1). 
For high temperatures T » Tk, as described elsewhere [25], the resistivity is well 
described by the Hamann result [36]. In the Fermi liquid regime T « Tk the inset in 
Figure [3] shows the expected behaviour. In the Kondo regime the values of the 
coefficient extracted from a least squares fit in the region < T < Q.ITk are 5.7, 5.8, 6.4 
and 6.6 in going from the symmetric case to = —2 A. These values agree to within 
8% of the exact result c = 6.088. Quantum Monte Carlo results for the symmetric 
case for [//vrA < 3 give the coefficient to within 19% of the exact result [H]. The 
current approach which is not restricted to very low or very high temperatures also gives 
accurate results for the resistivity in the cross-over region T ~ Tk [25]. The T = 
resistivity p(0) satisfies p(0) ~ sin^(7mrf/2) in accordance with the Friedel sum rule for 
the single-particle spectral density [26] . 

A second Fermi liquid relation relates the linear coefficient of the thermopower to 
the local level occupancy, ndiT = 0), and the linear coefficient of specific heat, 7, [29|[30] 

lim|^}=7rcot(7rM0)/2) (10) 

This is a much more difficult Fermi liquid relation to satisfy, because the local level 
occupancy rid^O) is obtained by integrating the single-particle spectral density up to the 
Fermi level, the linear coefficient of specific heat is obtained as a second derivative 
— lim I ^ Qr^P I of the Free energy and lim t^o | | involves evaluating the 
transport integrals using the finite T spectral densities pd{oJ,T). All these quantities 
are evaluated using the numerical renormalization group method and the above relation 
provides a very severe test of the accuracy of the method. 

The thermopower is shown in Figure H] for the mixed valent, empty orbital and 
Kondo regimes. The sign of the thermopower at low temperatures is determined by 
the sign of ( ) and is always positive for the range of parameters considered. 

At higher temperatures the behaviour of S{T) is more complicated. However, in the 
mixed valent and empty orbital regimes the single resonance always lies above the 
Fermi level (see Figure [1]), and from equation (7) the thermopower is positive at all 
temperatures. In the Kondo regime a distinctive maximum appears in the thermopower 
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Figure 1. The spectral density pd{uJ,T) in tlie mixed valent regime e^/A ~ —1 for 
several values of the reduced temperature t = T/Tk and for [//ttA = 4. Note that A 
is a more appropriate energy scale in this regime than Tk- For the present parameters 
Tk = 0.59A 



dX T Tk- At higher temperatures, the thermopower can become negative and for 
sufficiently high temperatures positive again. This bahaviour is associated with the 
temperature dependence of the sateUite peaks at and ed + U (see Figure [2]). The 
low temperature behaviour of the thermopower (inset to Figure H]) shows the expected 
linear in T Fermi liquid behaviour and one can use this to extract lim t^o | | • By 
calculating the linear coefficient of specific heat from the Free energy and ^(^(T = 0) 
from the spectral densities we can test the Fermi liquid relation (10). This is shown in 
Figure El 
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Figure 2. The spectral density pd{ijJ,T) in tlie Kondo regime e^/A — —3 for several 
values of the reduced temperature t — T/Tk and for [//vrA = 4. 



It is seen that limT-i-o li^s within 6% of cot(7rnrf/2) in nearly all cases. This 

agreement is remarkable and shows clearly that the method is capable of giving close 
to exact results for transport coefficients and finite temperature djTiamic properties in 
addition to the T = dynamics and the thermodynamics. We have not made any of 
the sophisticated corrections [31] to improve on the thermodynamics and only the even 
energy shell spectral densities have been utilized, hence there is scope for improvement. 
The characteristic peak aX T k, Tk in the thermopower in the Kondo regime is often 
referred to as the giant thermopower. This peak can become very large, although to see 
this one needs to keep Ud fixed and approach the Kondo regime by increasing U. In our 
calculations td was varied to approach the symmetric case for which the thermopower is 
zero as can be seen by setting = 1 in equation (10). From equation (10) it is clear that 
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Figure 3. The scaled resistivity p{T)/p{0) over the full temperature range in various 
regimes and for U/nA = 4, e^/A = 1 (□), e^/A = (A), e^/A = -1 (o), ea/A = -2 
(★), e^/A = —3 (+), Cd/A = —4 (x) and the symmetric case (•). It should be noted 
that a more appropriate energy scale in the mixed valent and empty orbital regimes is 
A. The inset is for (1 - (p(T)/p(0))) versus {T/TkY and shows the Fermi hquid 
behaviour for T « Tk in the Kondo regime. 



keeping fixed and increasing U leads to an exponential increase in the thermopower 
since 7 ~ l/T^- ~ e"^, a > thus giving a giant thermopower in the Kondo limit. The 
vanishing of the thermopower for half-filling is an artifact of the simplified model used 
here which neglects the non-resonant scattering of conduction electrons as discussed 
in [32]. An extension of the present calculations to include non-resonant scattering, 
leads to enhanced thermopowers for half-filled systems and describes the qualitative 
behaviour of the thermopower of concentrated Kondo compounds [33] . 
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Figure 4. The Thermopower S{T) over the full temperature range in various regimes 
and for U/nA = 4 e^/A = +1 (□), e^/A = (A), e^/A = -1 (o), e^/A = -2 (★), 
Cd/A = -3 (+), ed/A = -4 (x). The inset is for S{T)Tk/T versus T/Tk and shows 
the linear in T Fermi liquid behaviour for T << Tk- 



3. Conclusions 

In this Letter we have presented new results for the thermopower and resistivity of the 
Anderson model over the full temperature range of interest and in several parameter 
regimes by using the numerical renormalization group approach. The Fermi liquid 
relation for the linear coefficient of the thermopower is satisfied with remarkable 
accuracy, as is that for the coefficient of the resistivity. The Hall coefficient, thermal 
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Figure 5. liniT^o |^^^} (dashed line with diamonds) and cot(7rnrf/2) (solid line) 
versus Ud for the various parameter sets ranging from Kondo to mixed valent and 
empty orbital behaviour. The two curves lie within 6% of each other in all cases, 
except for the empty orbital case e^/A = +1 where the error is 10%. 



conductivity and finite temperature dynamic susceptibilities have also been calculated 
and will be presented elsewhere |26l [35] . Our results are in contrast to NCA calculations 
[121 [E] which violate Fermi liquid relations at low energies, and to perturbation in U 
calculations which although accurate eventually break down for U/ttA > 2.5 [12]. We 
are able to treat a much wider range of regimes than the Monte Carlo approach [TS] 
which has so far only been applied to the symmetric case where the thermopower is 
zero and which often builds in the sum rules which constitute independent tests within 
our approach. The techniques presented here for extracting the spectral densities of the 
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Anderson model provide a new highly accurate approach to the d — oo Hubbard model. 
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